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Abstract 

We study the stratum in the set of all quadratic differential systems x = 
P2{x,y), ij = Q2{x,y) with a center, known as the codimension-four case (54- 
It has a center and a node and a rational first integral. The limit cycles under 
small quadratic perturbations in the system are determined by the zeros of the 
first Poincare-Pontryagin-Melnikov integral /. We show that the orbits of the 
unperturbed system are elliptic curves, and / is a complete elliptic integral. 
Then using Picard-Fuchs equations and the Petrov's method (based on the 
argument principle), we set an upper bound of eight for the number of limit 
cycles produced from the period annulus around the center. 

Introduction 



The conditions for a plane polynomial quadratic vector field to have a center are 
known since the beginning of the last century (Dulac (1908), Kapteyn (1912)). In 
the space of all quadratic systems, the systems with a center form a union of four 
irreducible affine algebraic sets : 

• Hamiltonian (Q^) 

• reversible (Qf) 

• generalized Lotka-Volterra (Qf^) 
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• codimension-four set {Q4) 

(the subscripts indicate the co-dimension of each algebraic subset), see Zolg-dek [11] 
for a modern proof. Strictly speaking, the above classification concerns only the cen- 
ters themselves. There are quadratic systems with two centers, one of them in 
and the other in n Q^^ , see [7], p. 148. In the present paper we are interested 
in the maximal number of limit cycles which can emerge from the periodic orbits of 
a quadratic system with a center, after a small quadratic perturbation. The corre- 
sponding number is the cyclicity of the related open period anmihis. A quadratic 
center is said to be generic, if it does not belong simultaneously to two of the above 
algebraic sets. The cyclicity of the period annulus of a generic center depends on 
the number of the zeros of the first Poincare-Pontryagin-Melnikov function, while in 
the general case a higher-order analysis is needed. The corresponding higher-order 
functions were determined in Iliev [7]. 

The cyclicity of the open annuli in the generic Lotka-Volterra case (Qf^) has been 
studied by Zol§,dek [11], and in the generic Hamiltonian case {Q^) by Horozov-Iliev [6] 
and Gavrilov [4] (for the codimension-one generic cases from the bifurcation diagram 
of Q2 see [1] and the references therein). Almost nothing is known about the generic 
reversible case (Q^) and nothing about the generic codimension-four one {Q4,). In 
the present paper we place an upper bound on the cyclicity of the (unique) period 
annulus in the generic codimension four case {Q^}. 

Theorem 1. The cyclicity of the open period annulus surrounding the center of any 
generic codimension-four plane quadratic system is less than or equal to eight. 

The conjectural exact upper bound, as it is well known, is three [11, 7]. To 
the end of this Introduction we briefly sketch our proof. A codimension-four plane 
quadratic system with a center has a well known rational flrst integral H. Using 
this one may check that the generic complexified orbit {{x,y) e : H{x,y) = h} 
of the system is an affine elliptic curve, and the corresponding Poincare-Pontryagin- 
Melnikov function I = I{h) is a complete elliptic integral on it. This remarkable fact 
(it seems to have gone unnoticed by the specialists) is the starting point of the paper. 
Still, a direct application of standard techniques like in [10, 9, 6, 3, 4] does not work. 
Indeed, the Poincare-Pontryagin-Melnikov function involves differential of the third 
kind with residues algebraic in the parameter h (and not only polynomial, as in the 
usual cases) . On the other hand, it turns out that / satisfies a Picard-Puchs equation 
of the form 

M2 0L2 0Li(/) =0 

in which Li = h-^ — 1 and L2, M2 are suitable second-order Picard-Fuchs differential 
operators. The role of Li and L2 in our approach is to annihilate the residues of 
the differentials involved in I. It is easy to see that / has at most as many zeros as 
G — Li(I) on the interval of consideration. We shall prove further that the solution 
space of L2 is a Chebyshev space (Proposition 14). This on its turn implies that if 
i? is a function with k zeros, then each solution G of the non-homogeneous equation 
L2{G) = R has at most k + 2 zeros (Proposition 2). As R = L2 o Li{I) = L2{G) is in 
the kernel of the second-order Picard-Fuchs operator M2, it remains to show that R 
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has at most 6 zeros which is achieved in a standard way by making use of the Petrov 
method [10], see Proposition 4. 



2 Remarks about the codimension-four case 

In complex coordinate z — x-\-iy, the system corresponding to the generic codimension- 
four case with a center placed at the origin becomes 

z ^ -iz + Az^ + 2\z\^ + , Q;eC\E, \a\ = 2. 

It is well known that the codimension-four case has a first integral of the form 



n(x,y) 



[ipix,y)f 



where (p and i/j are polynomials of degree tree and two, respectively. Taking a — b+ic 
and denoting for short 

Y^cx-{2 + b)y, K = ^ > 1, 

one can easily find the explicit form of the first integral, namely 

[82/(1 + 1(1 + «:r3)]2 



[l-8y + kY^Y 



Therefore, ip = defines a trident curve and = is a parabola. Since (/^(O, 0) = — |, 
■0(0,0) = 1, the period annulus around the center is placed inside the domain fl 
determined by < < V- Taking ^ ijj ^ 1 - 8y + kY^, X > 0, then in 
(x, y) — {X, Y) coordinates, the system has in Q a first integral of the form 

H(x, y) = ^^^—^ l^-^i^f + i^f + (1 - x^)y - + 3 j ■ (1) 

(Compare to [7], where an additional rotation of the very initial coordinates takes 
place) . Although the change of the variables 

(x, y) ^ (x, y) 

does not define a bi-rational projective transformation (but a double ramified cover- 
ing), it induces a bi-rational map 

It is seen that the level sets of the first integral H{x,y) are cubic plane curves and 
hence they are (generically) eUiptic. 
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For convenience, from now on, the bars over the variables x, y will be omitted. 
The generating function I{t), whose zeroes correspond to limit cycles in the perturbed 
system, are given by the following complete elliptic integral (cf. [7], Theorem 2 (iii)) 

I{t)= / / x~^[ni + + l^3y^ + l^iiK^y^ - x^)]dxdy. (2) 

J JH{x,y)<t 

Our main purpose in this paper is to study how many zeroes the integral I{t) can 
have in the open interval corresponding to the period annulus around (1, 0). 



3 Picard-Fuchs systems 

In this section we derive several equations satisfied by the entries of (2). We also 
apply several reductions to express these integrals in the form we need. In particular, 
we calculate explicitly G, L2 and it! we mentioned in the introduction. 

We begin by introducing new variables y — yi — 1, h — 8(2 — b)t. Then the 
equation H{x, y) —t where H is given by (1) becomes (we will omit the subscript in 

H{x, y, h) = '^y^ - x^y - hx' - {k - l)y +'^{k - 1) ^ 0. (3) 

By using (3), it is not hard to verify that the period annulus around the center at 
(1, 1) corresponds to the levels h from the interval (— |, 
Let us denote for i, j e Z 



kAh) = [[ 

J JH{x,y, 



x^y^dxdy, h e 

h)<0 



3' S^/K J 



Then multiplying (3) by x^y^^^dx, respectively by x^^^y^dy and integrating, we obtain 
the identities 

|(j + 4)7,,,+3 - U + 2)/i+2j+i - h{j + l)7,+3j (4) 
-{k - l)(i + 2)7,,,.+i + ^{k - l)(j + = 0, 

|(i + l)/,,,+3 - (^ + 3)/,+2j+i - h{t + 4)/,+3,i (5) 

-(« - 1)(^ + l)J,,, + i + - 1)(^ + = 0. 

Multiplying (4) by i + 4 and (5) by j + 1 and subtracting, we come to 

K{i + j + + j + 5)/i+2j+i - (« - l)(i + 3j + 7)7i,,+i + 2{k - + = 0. 

In particular, for i = —6 and j = I one obtains /_6,2 = I-6,i- Therefore, the function 
in (2) takes the form 

I{h) = llj-efi + H2l-&,1 + A«3i-6,3 + - /-2,o) (6) 
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with all constants /Xj e M independent. 

Let us apply to (3) and (6) the change of variables {x, y) yx~^). Then (3) 

reduces to 

H{x, y) = -{K- l)a;3 - {k - l)x^y + ^y'' - y = h, (7) 
Ii,j{h) becomes —I-i-j-sj{h) and (6) becomes 

Hh) =11 (l^ix^ + l^2X^y + l^3y^ + 1^4 ( ^^ ) dxdy. (8) 

J JH{x,y)<h V V ^ / / 

By (7), H{x,y) = -~H{—x, —y), therefore the phase portrait of the related Hamilto- 
nian system has a central symmetry with respect to the origin. 
Next, we can use the following identities: 

r _ r _ 3^ r t 

h,2 — h,i — '° ~'~ '° ^ 5 

^ 10(«:-l)^°'° + + 5(^^^°'^' 
_3{K+l)h , + 6 

- 10k + + 
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to transform (8) into 

I{h) = fiihlo^oih) + fi2h,o{h) + f^3h,i{h) + /i4[2/-i,o(/i) + 3«;/i7_i,i(/i)]. (9) 

Following the standard way [6], one can derive a Picard-Fuchs system for the entries 
in (9). Its explicit form is as follows. 

i0,Q - ^-'0,0 + -'0,1' 

7-2/' +hi' I 

r -^r 4--T' ^-T' r 

-'1.1" g -'0,0 2 1.0 4 0.1 4 1.1' 
7_i,o = 3/i7ii,o + 2/ii,i, 

^-1.1 - ~;^^i.o + ~^-i,o + Y^-i.i- 
By using the above system, we see that 

hi' -I^ G{h), G{h) = (//i/i' + //3)/oo + /^2/; 1 + //4[-4/ir 10 + - 4)^1 1]. 
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Therefore 



I{h)^hJ^J-'G{OdC 



and I{h) has at most as much zeroes as G{h) in (— |, —■^^)- By the same system, 
the integrals I'q q and I[ i satisfy 

-3^/1/^,0 = (9«/i'-4)/^,o-4(«:-l)/r,i, .^q. 
-3KhI[^, = {9Kh'-4){I'^^,-I'l,), 

and the integrals I'_i q and i satisfy 

T' _ „ ^ II 
-'-1,0 — 2 -1.0 -1.1' 

^-1,1- ^^-1,0 ^1,1- 
Hence, the integral J = — 4/i7ii q + (Sk/i^ — 4)7^1^^ satisfies the second-order equation 

L2(/i) J = ^(k - \)\h{%Kh^ - 4)7('^ + (Gk/i^ + 8)7('_J, 

with 

Uih) = 5Kh - {9Kh^ - 8) — + h{9Kh^ - 4) 



dh ' 

Therefore, a similar equation L2{h)G — R (with a right-hand side R depending 
linearly on fii, 7q q and I[ i ) also holds. To calculate R explicitly, we first use (10) to 
obtain the identities 

(9/i2-4)(9«/i2_4) 

9h^-A 

III _ 32AKh^ + {72k - 108)/i2 - 48 12(k - 1)[243k,/?.4 - 36(k, + l)/^^ - 16] 

~ (9/l2-4)2(9/t/l2-4) (9/l2 - 4)2(9/^/i2 _ 4)2 1,1' 

_ 27/?.2 + 12 , _ lQ2Kh* + (144k, - 108)/?.^ - 48 , 
^'^ ~ (9/i2 - 4)2 (9/i2 - 4)2(9k/i2 - 4) 

A direct calculation then yields 

h[{ao + ai/i' + 02/^' + a3/^^)/o,o + (&o + hh'' + &2/^')/l,i] 

(9/i2 - 4)2(9«/i2 _ 4) ' ^^^^ 

with some constants Oj, 6^ depending linearly on ^i. Below, we shall use the explicit 
formulas for G, 7?, L2 just derived in order to prove our main result. 
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4 Proof of Theorem 1 

The proof of Theorem 1 follows from the next four Propositions, the first two of 
them being probably known. Let y be a finite-dimensional vector space of functions, 
real- analytic on an open interval (a, 6). 

Definition 1. We say that V is a Chebyshev space, provided that each non-zero 
function in V has at most dim{V) — 1 zeros, counted with multiplicity. 

Let S be the solutions space of a second-order linear analytic differential equation 

x" + ai{t)x' + a2{t)x ^ (12) 

on an open interval (a, b). 

Proposition 1. The solution space S of (12) is a Chebyshev space on the interval 
{a,b) if and only if there exists a nowhere vanishing solution Xq E S {xQ{t) ^ 0, 
\/t e (a, 6)). 

Remark 1. The question of existence of a non- vanishing solution is a recurrent 
question in many papers concerning zeros of Abelian integrals, see e.g. [9, 5, 2]. A 
Chebyshev space V in our sense is sometimes called an extended Chebyshev space, and 
it is said to be an extended complete Chebyshev space, provided that it has a complete 
flag of extended Chebyshev sub-spaces, see e.g. [8]. In the case when dim(\/) = 2 
the Chebyshev space V (in our sense) is an extended complete Chebyshev one if and 
only if it has a nowhere vanishing function. Therefore the notions of Chebyshev space 
(in our sense) and extended complete Chebyshev space (in the sense of [8]), as far as 
applied to the solution space of (12), coincide. 

Proposition 2. Suppose the solution space of the homogeneous equation (12) is a 
Chebyshev space and let R be an analytic function on (a, b) having k zeros ( counted 
with multiplicity). Then every solution x{t) of the non-homogeneous equation 

x" + ai{t)x' + a2{t)x = R{t) (13) 

has at most k -\-2 zeros on (a, b). 

Proposition 3. The solution space S associated to the differential operator 

d 

L2(/i) =5k/i-(9k/i2-8)— + /i(9k/i2-4) — , k>1 (14) 



is a Chebyshev system on the interval (— oo, — 




Let R — L2 o Li{I), where / is the Abelian integral (9) and Li — h-^ — 1. A 
suitable for our purposes expression for R is obtained in (11). 

Proposition 4. The Abelian integral R{h) has at most 6 zeros (counted with multi- 
plicity) on the interval ( — |, —3^)7 k,> 1. 
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Proof of Theorem 1, assuming Propositions 1-4. The Abelian integral Li{I), 
I given by (9), is a solution of the non- homogeneous equation L2{G) = R. According 
to Proposition 4, Proposition 3 and Proposition 2 the integral Li[I){h) has at most 
8 zeros on the interval (— oo, —3^)- The integral I{h) has the same number of zeros 
as Li{I){h) on the same interval. Finally, the functions (9), (8), (6) and (2) have the 
same number of zeros in the respective intervals. □ 

To the end of the paper we prove the above Propositions 1-4. 



Proof of Proposition 1. Let (xi, X2) be a fundamental set of solutions of (12) and 
consider the map 

p : (a, b)^ ^ P^M : t ^ [xi{t) : X2it)]. 

As the Wronskian of Xi , X2 is non- vanishing, then the map p is non-degenerate {dp{t) ^ 
0) and hence monotonous. 

The solution space of (12) is Chebyshev if and only if the map p is injective. The 
solution space of (12) contains a nowhere vanishing solution if and only if the map p 
is not surjective. 

As the circle is not homeomorphic to an open interval, then the monotonous dif- 
ferentiable map p cannot be surjective and injective at the same time. It follows that 
if p is injective then it is not surjective. If, on the contrary, p is not surjective, then 
the monotonicity of p implies that the image of (a, b) under p is an open subinterval 
of and p is injective. □ 



Proof of Proposition 2. Let (xi, X2) be a fundamental set of solutions of (12), such 
that xi{t) is a nowhere vanishing solution. The change of the variables x — > x/x\{t) 
does not change the number of the zeros of the solutions of (12), which is transformed 
to a linear equation with a fundamental system of solutions {1, ^^}- As the vector 

space spanned by Xi, ^2 is Chebyshev, then the function is strictly monotonous 

on (a, 6). The change of the independent variable t — > r = is therefore regular 
and transforms further the above linear equation to an equation with a fundamental 
system of solutions {1,t}. Therefore the corresponding differential operator is a 
multiple of More precisely, the regular change of variables 

X X2{t) 

(x, t) ^{y,T), y = — -, T = — - 

Xi{t) Xi{t) 

transforms equation (13) to 

and hence each solution of the non- homogeneous equation (13) has at most k + 2 
zeros on (a, 6) (counted with multiplicity). □ 
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Proof of Proposition 3. Let {S{h) : h e (— |, —3^)} be the continuous family of 
periodic orbits defined by {H = h}, with H in the form (7). Then G{h) = 



where u; is a linear combination of elliptic differentials of the first and second kind. 
Therefore the residues of uj are solutions of L2. The only residues of ut are at (0, 
where y is one of the roots of — y — h and they are easily computed: 

-4h + {3Kh^ - 4)y 

Res{Q^y)Uj = — . 

^ Ky^ — 1 

For h < the polynomial ||/^ — y — h has one real root which we denote by y^. 

We shall show that the solution f{h) = Res(o,yo)^ of L2X = does not vanish in the 
interval (—00, — ^^). Indeed, on this interval yo{h) is a strictly increasing function 

and yo{h) < —\J\ = l/o(~3^)- It remains to show that —Ah + {Suh'^ — 4)7/0 7^ 0. 
The identity 

-Ah + (SKh^ - 4)yo = -4(|yo' - ^o) + (3«/i^ - 4)yo = Kh{3h^ - ^yl) 
implies that f{h) = on {—00, —3^) if and only if /i = 2yo/3. Now 

gives yo — '^\J\- Vo — yo{h) is a strictly increasing function in (—00, —3:^] 

and yo(— 3:^) = "^/f which is the needed contradiction. Thus the solution space of 

L on (—00, —3^) contains a nowhere vanishing function and hence is a Chebyshev 
system. □ 

The above result cannot be improved, as shown by the example L — 

Proof of Proposition 4. According to (11), it suffices to show that any hnear com- 
bination of the form Ps{h^)lQ Q{h) + Q2{h'^)I'i^i{h) where P3, Q2 are real polynomials of 
degree at most three and two, has at most 6 zeros. We note that Iq q, ![ 1 are complete 
elliptic integrals of the first and second kind respectively, satisfying the second-order 
Picard-Fuchs system (10). 

We introduce a new variable s e (1, «;) through h = —\\/s/k and denote by dot 
the differentiation with respect to s. Also, denote for a convenience Ji(s) = lQQ{h{s)), 
J2(s) = I[ -^{h{s)). By (14) and (11), we obtain the equation 



LoG 



_1 d _ 5' 
■'d^~2d~s~36 



P3{s)Jl{s) + Q2{s)J2{s) 
{S - k)2(s - 1) 



Hence, we will need information about the zeroes of the right-hand side in the interval 
(1, k) C (1, 00). Equation (10) implies that J(s) = {Ji{s), J2(s))^ satisfies the system 
of hypergeometric type 
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or equivalently 



6{s-1){s-k)J{s)= {] ' " ] \J{s) (15) 

1 — s s — 1 



Let us consider for any n eN the vector space 

K = {Pn-Jl + Qn-l-Ji ■■ Pn, Qn-1 £ ^[^1, deg P^, Qm < m}. 

Proposition 4 follows from the following more general result. 

Proposition 5. The vector space Vn is Chebyshev on the interval (1, k) : each element 
has at most dimVn — 1 = 2n zeros (counted with multiplicity). 

Proof of Proposition 5. We use the Petrov method in the complex domain V = 
C \ (— oo, 1), see [10, 9]. The characteristic exponents of (15) at 1, k, oo are equal to 
{0, 0}, {0, 0}, { — |, |}, respectively. The function Pn-Ji + Qn-i-J2 is holomorphic in 
a neighborhood oi s — k, (this value corresponds to the center of the system dH — 0, 
with H the symmetric Hamiltonian given by (7)), and has a logarithmic singularity 
in a neighborhood of s = 1 (which corresponds to the saddle point of the symmetric 
Hamiltonian system) . The function Ji is a complete elliptic integral of the first kind 
and therefore does not vanish. Consider the function 

Ms) 

which is real-analytic in the complex domain T>. We apply the argument principle to 
the smaller domain 

V^ = Vn{s:\s-l\>e}n{s:\s\< 

For this purpose, we consider the increase (or decrease) of the argument of F when s 
makes one turn along the boundary of in a positive direction. The following facts 
are easily deduced from the asymptotic expansions of J near the singular points of 
the Fuchs system (15). 

1. Along the boundary of the small disc {|s — 1| — s} the increase of the argument 
of F is bounded by a value close to zero. 

2. Along the boundary of the big disc = ^} the increase of the argument of F 
is bounded by a value close to 27rn = 27rmax{n, n — 1 + |}. 

3. Along the interval (— oo, 1), the imaginary part of F equals 

^ /XT M-s) ^ , ,detW(s) 

where 

' Ji(s) Ji(s) 



W{s) -- 

is a fundamental matrix of (15). 



J2{S) J2{S) 
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4. The determinant of the fundamental matrix 1^ is a rational function in s and 
in fact a non-zero constant. 

Summing up the above facts wc conclude that the increase of the argument of F 
along the boundary of is bounded by 2n. Therefore F, and hence P„(s). Ji(s) + 
Qn-i{s).J2{s) has at most 2n zeros (counted with multiplicity) in V, and hence in 
(!,«;). Proposition 5, and hence Proposition 4 are proved. □ 
This also finishes the proof of Theorem 1. 
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